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Abstract. We present baseline calculations of initial-state shadowing and final- 
state absorption effects on J /if) production in nucleus-nucleus collisions at the 
Relativistic Heavy Ion Collider. We show predictions for Au+Au and Cu+Cu 
collisions at \JSnn — 200 GeV and Cu+Cu collisions at ^/Snn = 62 GeV as a 
function of the rapidity, y, and the number of binary nucleon-nucleon collisions, 
AUi- 
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The nuclear dependence of J r /ip production has been studied extensively in 
fixed target nuclear collisions [H] where cold matter effects were attributed solely to 
nuclear absorption. However, the nuclear parton distributions are modified relative 
to those of free protons. At V~S = 17.3 GeV, the modification, referred to here as 
shadowing, gives a modest enhancement of J/ip production at midrapidity since the 
J /ip is produced in the antishadowing region [El- As Vs increases, the momentum 
fractions, x, probed by J/ij) production at y = decrease, leading to stronger 
shadowing effects at collider energies. New measurements at the Relativistic Heavy 
Ion Collider (RHIC) on nuclear collisions may determine the importance of dense 
matter effects on «//•</> production. For these analyses to be meaningful, it is essential 
to have a proper baseline for quarkonium suppression in AA collisions due to cold 
matter effects. 

This paper studies the interplay of shadowing and absorption in nuclear colli- 
sions at RHIC. We address both spatially homogeneous (minimum bias) results as 
a function of rapidity and inhomogencous (fixed impact parameter) analyses as a 
function of the number of nucleon-nucleon collisions. Our calculations, described 
in detail in Ref. [0] for d+Au collisions, agree relatively well with the preliminary 
PHENIX d+Au data [H]. Thus if J/ip production is unaffected by quark-gluon 
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plasma formation, these predictions should describe the rapidity and centrality de- 
pendence of J/tp production in A A collisions. Since RHIC has completed both 
Au+Au and Cu+Cu runs at ^S NN = 200 GeV and a ^S NN = 62 GeV Cu+Cu 
run, we present predictions for shadowing and absorption effects on J /if) production 
in these three systems. 

Our J/ip calculations generally employ the color evaporation model (CEM) [ 
E|- The J ftp rapidity distribution in an AB collision at impact parameter b is 

xSt s (r,z)Si s (b-r,z>) . (1) 

where — (M/ \J Snn) exp(±y) at leading order. The parameter Fm is the 
fraction of cc pairs below the DD threshold that become J/ip's [E]- Since the 
shadowing ratios are independent of the order of the calculation [0 [7], the AA/pp 
ratios are calculated at leading order to speed the numerics. We use m c = 1.2 GeV 
and Q = M [El with the MRST LO parton densities [Ej. 

The parton densities in the nucleus, Ff-(x, Q 2 ,r, z) = pa{s)S p S (A, x, Q 2 7 r, z) 
ff(x,Q 2 ), are the product of the nucleon density [Gj], pa(s), the nucleon parton 
density, f^(x,Q 2 ), and a shadowing ratio, S P S (A, x, Q 2 , f, z), where r and z are 

the transverse and longitudinal position of the parton with s = \Jr 2 + z 2 . The two 
subscripts on Sp g refer to the shadowing parameterization and the spatial depen- 
dence, respectively. Most of our results are presented for the EKS98 shadowing 
parameterization [^3 but we also compare to predictions with the three FGS pa- 
rameterizations [lll|: FGSo, FGSh and FGS1. Since A = 63 is not available from 
FGS, the results here are calculated with A = 40. 

Figure ^a) compares the four homogeneous gluon ratios, 5g KS and 5*p GSi for 
Q = 2m c . The thin curves are for A = 197 while the thick curves show A = 63 
(A = 40 for FGS). FGS predicts more shadowing at small x. The decrease of 
shadowing with A is stronger for FGS since a smaller A is used. The EKS98 
and FGS1 ratios for the lighter A are very similar for 0.001 < x < 0.2. While 
all parameterizations show significant antishadowing, the EKS98 antishadowing x 
range is larger. Shadowing alone will give an effective A dependence as a function 
of rapidity with y > corresponding to low X2, effectively mirroring the curves in 
Fig. ^a) for dA while, in AA collisions, the result should be similar to the product 
of the curves with their mirror images. 

Figure H^b) shows the average value of X2 for ^/Snn = 200 (solid) and 62 
(dashed) GeV. At midrapidity at 200 GeV, (xz) ~ 0.01, near the point where 
Sp < 1. In the forward region, (X2) ~ 10~ 3 at y ~ 2, clearly in the effective low x 
regime while for negative rapidity, y ~ —2, (#2) ~ 0.1, in the antishadowing region. 
At 62 GeV, the average (2:2) is larger and the rapidity distribution is narrower. 
At midrapidity, (x 2 ) ~ 0.05, in the antishadowing region. Thus, the predicted 
shadowing effects should be quite different for the two energies, as we will show. 
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Fig. 1. (a) The gluon shadowing parameterizations at scale fi = 2m — 2.4 GeV for 
EKS98 (solid), FGSo (dashed), FGSh (dot-dashed) and FGS1 (dotted). The thin 
curves are for Au while the thick curves are for Cu (Ca for the FGS parameteriza- 
tions) . (b) The average value of X2 in pp collisions as a function of y for ^fS = 200 
(solid) and 62 (dashed) GeV. 

We assume that the inhomogeneous EKS98 and FGSo parameterizations are 
proportional to the parton path through the nucleus, Sp normalized so that 

(1/A) J d 2 rdzp A {s)S Vp {A,x,Q 2 ,r,z) = S l P {A,x,Q 2 ). The inhomogeneous FGS 
parameterizations are normalized so that {1/A) J d 2 sTA(s)S FGShl (A,x,Q 2 ,s) = 

^FGSh,i(^' x i Q 2 )- 

The survival probability, S , ^| bs , for J/ip absorption by nucleons is S abs (r,z) — 
exp{— J dz' pA(r, z')a a b s (z' — z)} where z is the production point and z' is the 
absorption point. We consider both color octet and singlet absorption. The octet 
\(cc)gg) state can convert to a singlet with a formation time of 0.25 fm at negative 
rapidity. Otherwise if the J/ip, 4>' and \c hadronize outside the nucleus, they have 
identical, finite octet cross sections. The singlet cross section grows quadratically 
with proper time until the formation time of each charmonium state. Thus the 
individual charmonium states have different asymptotic cross sections depending 
on their relative radii but absorption is ineffective if the state forms outside the 
target. See Ref. [E| for more details. 

In Fig. [21 we contrast the AA/pp ratios in Au+Au and Cu+Cu collisions cal- 
culated with EKS98 and asymptotic absorption cross sections, a ahs — 0, 1, 3 and 
5 mb, for octet absorption (upper plots) and color singlet absorption (lower plots). 
The zero value illustrates the effects of shadowing alone. A 3 mb absorption cross 
section was needed to obtain agreement with the E866 800 GeV J/ip data [EI at 
xf ~ [E]- The absorption cross section has been predicted to both decrease [ 
IT5] and increase [EB] with energy. The PHENIX d+Au data [|4j suggest er abs < 3 
mb at RHIC [EjJ, in accord with a constant or decreasing <7 a bs- 

The AA/pp ratios are symmetric around y = 0. The effect of octet to singlet 
conversion at large y is negligible for \?Snn = 200 GeV when cr a b s > 2.5 mb, the 
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singlet cross section assumed in the model [E|. The results with singlet absorption 
differ from those with shadowing alone only at \y\ > 1—2 at 200 GeV and \y\ > 0.5—1 
at 62 GeV. The larger values of cr a b s correspond to the greater rapidity range. 
Formation is within the 'target' at negative rapidity and within the 'projectile' at 
positive rapidity. 
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Fig. 2. The AA/pp ratio with the EKS98 parameterization as a function of y for 
octet (upper) and singlet (lower) absorption. In (a) and (b) we show the Au+Au 
results at 200 GeV while the Cu+Cu results are shown at 200 GeV (c) and (d) 
as well as at 62 GeV (e) and (f). The curves are a a bs = (solid), 1 (dashed), 3 
(dot-dashed) and 5 mb (dotted). 

There are two peaks in the 200 GeV ratios at y w ±1.5, the location of the 
antishadowing peak at x ps 0.1. There is a dip at y = where the shadowing region 
in x of one nucleus coincides with the EMC x region of the other. The convolution 
of the two shadowing ratios causes AA/pp < 1 over all y rather than a ratio greater 
than unity in some regions, as in dAu/pp [|2|. The results are similar for Au+Au and 
Cu+Cu except for the relative magnitude since shadowing and absorption effects 
are both reduced for the smaller A. 

On the other hand, at \/Snn — 62 GeV, the antishadowing regions of the two 
nuclei coincide, giving a peak with AA/pp > 1 for y w 0. Away from midrapidity, 
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the ratio decreases until y ~ ±2.5 where x is in the 'Fermi motion' regime, ap- 
proaching the edge of phase space. Note also that at the lower energy the J/tp can 
be absorbed inside the target over a larger range of y since the effective formation 
time is shorter at the lower energy. 

Figure [3] compares the shadowing parameterizations for the three systems as- 
suming octet absorption with cr a b s = 3 mb. The differences are largest for Au+Au, 
the only system where the EKS98 and FGS parameterizations can directly be com- 
pared since A = 197 is an option for them all (except FGSo which takes A = 206). 
The differences between the shadowing parameterizations are larger in AA than in 
dA collisions since the combination of two nuclei enhances the differences. Since 
the parameterizations are most similar in the antishadowing region, the range of 
predictions is reduced at V<Saw = 62 GeV. 
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Fig. 3. The ratio AA/pp as a function of y for octet absorption with <7 a bs = 3 
mb and the EKS98 (solid), FGSo (dashed), FGSh (dot-dashed) and FGS1 (dotted) 
parameterizations for Au+Au at 200 GeV (a) and Cu+Cu at 200 GeV (b) and 62 
GeV (c). 

Finally, the spatial dependence of the AA/pp ratios are shown in Fig. Q] for 
the three systems as a function of the number of binary nucleon-nucleon collisions, 
iVcoii, assuming cr a b s = 3 mb. We have chosen two rapidities to illustrate the results: 
y = (dashed) and 2 (solid), corresponding to the central and forward regions 
covered by PHENIX lepton masurements. We do not show the equivalent backward 
region, y = — 2, since they are identical to those at y = 2. The EKS98 results 
shown here are very similar to the FGS predictions, even those with inhomogeneous 
parameterizations. Note that if a smaller <7 a bs is required by the rapidity-dependent 
data, the ratios shown here would be closer to unity. One might naively expect that 
the ratios at y = 2 would be lower than those at y = since xi is lower than at 
midrapidity but because y = 2 is in the antishadowing region of x\ at 200 GeV, its 
ratio is higher. On the other hand, at 62 GeV and y = 2, X2 is in the shadowing 
regime while x\ is in the EMC regime, making the ratio lower. In the backward 
region, x-i and X\ are interchanged but the results remain identical. 
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Fig. 4. The ratio AA/pp as a function of iV co ri for a 3 mb octet absorption cross 
section and the EKS98 parameterization at y = (dashed) and y = 2 (solid) for 
Au+Au at 200 GeV (a) and Cu+Cu at 200 GeV (b) and 62 GeV (c). 

The results are shown as a function of N co n since J/tp production is presumed 
to be a hard process. We have not shown the ratios for b > 2.\Ra where iV co ii ~~ * 1- 
The number of nucleon-nucleon collisions in Au+Au is considerably larger than the 
maximum iV co n for Cu+Cu. Since the inelastic nucleon-nucleon cross section is 
lower at 62 GeV, iV co ii is somewhat reduced relative to 200 GeV. The results are 
similar as a function of the number of participating nucleons, N pai t , the centrality 
variable for soft particle production. We prefer to discuss J/tp production as a 
function of N co \\ or N palt since these variables depend only on the nuclear density 
distributions. No additional model assumptions about the nature of the medium 
are needed, as is the case for other variables such as the path length, L, [^] or the 
energy density, e [ITU], 

PHENIX has shown that these cold matter effects are important in d+Au col- 
lisions [ 4 . If there are no additional effects on J/tp production in AA collisions, 
the predictions shown here should describe the AA results. If quark-gluon plasma 
production causes further J/tp suppression in AA interactions, then there should be 
deviations from this pattern. The extent of these deviations would depend on the 
rapidity range of the produced plasma. Since dense matter effects seem important 
for other ^4^4 observables some additional J /if) suppression might be expected. 
However, recent calculations of lattice-based potential models [ I21| and J/tp spec- 
tral functions [ [221 show that the J/ip itself might not break up until T ~ 2.5T C . 
On the other hand, the -0' and \c contributions, 40% of the yield, should be re- 
moved if these states break up at 0.7 < T/T c < 1.1, as expected. The loss of these 
contributions should be observable in the rapidity dependent ratios as additional 
absorption, especially at midrapidity, relative to the dAu/pp ratios shown in Ref. [ 
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